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Abstract 

Hadronic production of the doubly charmed baryon H cc and H+) is investigated under 

the general-mass variable-flavor-number (GM-VFN) scheme. The gluon-gluon fusion mechanism 
and the intrinsic charm mechanisms, i.e. via the sub-processes g + g — > (cc)[ 3 <Si]g + c + c, 
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\Q ■ g + 9 -> (cc)^ So} e + c + c; g + c^ (cc)[ 3 Si]3 + c, g + c -> (cc)^S ] 6 + c and c + c -> (cc) [ 3 5i]g + g, 

O ■ 

c + c — > (cc)[ 1 5o]6 + 5, are taken into account in the investigation, where (cc)[ 3 .Sijg (in color 3) and 
(cc)f 1 So]q (in color 6) are two possible iS-wave configurations of the doubly charmed diquark pair 
(cc) inside the baryon S cc . Numerical results for the production at hadornic colliders LHC and 
TEVATRON show that both the contributions from the doubly charmed diquark pairs (cc)[ 1 S'o]6 

•i-H . 

and (cc)[ 3 /Si]g are sizable with the assumption that the two NRQCD matrix elements are equal, 
and the total contributions from the 'intrinsic' charm mechanisms are bigger than those of the 
gluon-gluon fusion mechanism. For the production in the region of small transverse-momentum 
Pt, the intrinsic mechanisms are dominant over the gluon-gluon fusion mechanism and they can 
raise the theoretical prediction of the H cc by almost one order. 
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I. INTRODUCTION 



The heavy hadron S+ may have been observed by SELEX Collaboration already [1, 2], 
although some comments [3] pointed out that the measured lifetime is much shorter and the 
production rate is much larger than most of the theoretical predictions [4-8]. It is predicted 
that at the fixed target experiment, only about 1CT 5 of A+ events in its total sample are 
producted by S+, however the SELEX collaboration has found that almost 20% of A+ events 
in its total sample are producted by S+ . 

In the literature, most of the perturbative QCD (pQCD) calculations and predictions 
for S cc 1 hadroproduction are based on the 'gluon-gluon fusion mechanism' i.e. that via the 
subprocess, g + g — > (cc)[ 3 Si]s + c+c only. Whereas, the subprocess g + g — > (cc)[ 1 £ , ]6 + c+c 
may also contribute to the production [9]. It is because that and S+ + contain the higher 
components {ccqg) (here q — u,d) etc in their Fock space expansion, so the corresponding 
subprocess should be taken into account. 

The discussion shown in Ref. [9] indicates that an inclusive production rate of or S+ + 
can be factorized into two parts, one part is to produce two free c quarks, which can be 
calculated by pQCD, another part is to make these two free c quarks into a cc diquark pair: 
(cc)g[ 3 S'i] or (cc^pSo], then the diquark pair hadronizing either into by absorbing a 
quark d or into S+ + by absorbing a quark u for (cc)g[ 3 S'i], or either into S+ by absorbing 
a quark d or into by absorbing a quark u but both absorbing an extra soft gluon for 
(cc) 6 [ 1 5'o], all of which can be attributed to non-relativistic QCD (NRQCD) matrix elements 
[10]. In most of the existent calculations for the hadronic production of S cc , the cc diquark 
pair is assumed to be in 3 S\ configuration and in the color representation 3 ((cc)^ 3 ^]). 
Whereas according to power counting in velocity v c , the velocity of the heavy c quarks in the 
baryon, the NRQCD matrix elements h\ and (defined in Eq.(3)) for the nonperturbative 
transition, which correspond to the two configurations of the diquark pair {hi is that for 
(cc)g[ 3 S'i] and h 3 is that for (cc) 6 [ 1 S , ]), are at the same order of v c [9]. Hence to give a full 
estimation of the hadronic production of S cc , we think that (cc) g^Si] and (cc^^So] should 
be treated on the equal footing. 

Moreover, as pointed out in Refs.fll, 12], the so-called 'intrinsic' charm mechanism can 

1 Throughout the paper, S cc denotes S+ or i.e., the isospin-breaking effects are ignorable here. 
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give sizable contribution to the charmonium hadroproductionfll], and to the B c hadropro- 
duction, especially in small p t region[12]. Therefore, in addition to considering two configu- 
rations of diquark pair in different color representation 3 and 6 for the gluon-gluon fusion 
mechanism, it is also interesting to see how important of the 'intrinsic' charm mechanisms 
via the sub-processes g + c — > (cc) [ 3 <Si]jj + c, g + c — > (cc) I 1 So] 6 + c and c + c — > (cc) pSiJg + g, 
c + c — > (cc)[ 1 S'o]6 + <7, in hadronic production of and H+ + precisely. 

Principally, the 'intrinsic' charm mechanism induced by the heavy charm quark is greatly 
suppressed by the parton distributions in comparison with the valance and the sea of light 
quarks and gluon, but it is 'compensated' by 'greater phase space' and lower order of in- 
teraction coupling of QCD. Namely the 'intrinsic' processes are 2 — > 2 sub-processes at the 
order of 0(a^), while for the gluon-gluon fusion subprocess, its leading contribution starts 
at 0(ctf) and is a 2 — > 3 process. 

For a fixed target experiment which can reach to the region of very small transverse mo- 
mentum p t , the production of the doubly charmed baryon S cc should additionally involve 
more 'mechanisms', such as the mechanisms of the so-called intrinsic charm fusion with the 
subprocesses c + c — > (cc) pSijg and c + c — > (cc)! 1 Sole, which contribute to the produc- 
tion only with very small p t but whose nature essentially is non-perturbative for QCD. The 
theoretical predictions on the hadronic production rate all can be based upon the perturba- 
tive QCD (pQCD) only, though the existent ones are orders of magnitude smaller than the 
SELEX observation as pointed out by Ref.[3]. Nevertheless, we think that it is worthwhile 
to consider more mechanisms than that in the existent predictions, and to use the updated 
parton distribution functions (PDFs) in the general-mass variable-flavor-number (GM-VFN) 
scheme to re-estimate the S cc hadroproduction so as to cover a so widen p t region as pQCD 
is applicable. Especially, more attention to the so-called intrinsic charm production mech- 
anism, that is through the subprocesses g + c — > (cc) pSijg + c, g + c — > (cc)[ 1 S'o]6 + c and 
c + c — > (cc)[ 3 5'i]3 + g, c + c — > (cc)! 1 Sole + g, should be payed. 2 

2 The reliable estimate of the production so far can be that in terms of pQCD only, so here we take into 
account all the mechanisms which are calculable by pQCD. Therefore, here the so-called intrinsic charm 
fusion with the subprocesses c + c — > (cc)[ Si] 3 and c + c — > (cc)[ So]e are n °t considered (because 
they are of non-perturbative QCD as mentioned above). Since the 'higher order' mechanisms with the 
subprocesses: c+c — > (cc)[ 3 Si]g+.g, c+c — > (cc)[ 1 So]6+3 are also taken into account so as to 'complete the 
estimate', so in order to guarantee pQCD applicable and the obtained results being reliable, we compute 
the production always to put on a sizable cut on the transverse momentum p t of the produced (cc)-pair. 
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This work is devoted to give a comparative studies of various production mechanisms, and 
is also served as a cross-check of the pQCD calculation for the gluon-gluon fusion mechanism, 
because the results given in Ref. [5] and Ref.[6, 7] are in disagreement. Our results satisfy 
the gauge invariance at the amplitude, and our results agree with that of Ref. [5] except for 
an overall factor 2. 

When combining the results of 'intrinsic' charm mechanism with the gluon-gluon fusion 
mechanism, one needs to make some subtractions to the 'intrinsic' mechanism so as to avoid 
'double counting'. To perform the subtraction, we adopt the general-mass variable-flavor- 
number (GM-VFN) scheme [13-15], in which the heavy-quark mass effects can be treated 
in a consistent way both for the hard scattering amplitude and the PDFs. Moreover, it 
will be necessary to use the dedicated PDFs with heavy-mass effects included, which are 
determined by global fitting utilizing massive hard-scattering cross-sections. For instance, 
for the present analysis, the up-dated one CTEQ6HQ [16] is used. 

In Ref. [9], the production of E cc at e + e~ collider is treated carefully and hadronic produc- 
tion is estimated roughly by comparing with c-quark jet both by taking the fragmentation 
approach. In the present paper, alternatively, we take the full pQCD approach to do the 
estimate with more mechanisms, because we think the fragmentation approach becomes 
reliable only at the high p t regions where the fragmentation mechanism is dominant and 
also the results from the fragmentation approach show a strong dependence on the input 
parameter values [9]. At last, our results show that when assuming hi = /13, the contribution 
to the Hadronic production of S cc from the doubly charmed diquark pair in (cc^^So] can 
be sizable as that from (cc)g[ 3 S'i]. 

The paper is organized as follows. In Sec. II, we shall first give the formulation for the 
hadronic production of H cc within the GM-VFN scheme, and then present in some more 
detail the formulae for both the gluon-gluon mechanism and the 'intrinsic' charm mechanism. 
In Sec. Ill, we present the results for the subprocess and make a comparison with those in 
the literature. In Sec. IV, we present the numerical results for the hadronic production of 
H cc and make some discussion over them. The final section is reserved for a summary. 
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II. FORMULATION UNDER THE GM-VFN SCHEME 



Under the general-mass variable-flavor-number (GM-VFN) scheme [13-15], according to 
pQCD factorization theorem the cross-section for the hadronic production of H cc can be 
formulated as below: 



a = Ffj i (x 1 ,/j 1 )Ff l2 (x 2 ,n)Q^a gg ^ Scc (x 1 ,X2,n) 

'&cc-+H cc (x 1 ,X 2 ,tl)-\ , (1) 



where the symbol • • • means even higher order a s terms. F H -(x,fj,) (with H = Hi or H 2 ; 
x = xi or x 2 ) is the distribution function of parton % in hadron H . a stands for the cross- 
section of the corresponding subprocess. For convenience, we have taken the renormalization 
scale /j,r for the subprocess and the factorization scale fj,p for factorizing the PDFs and the 
hard subprocess to be the same, i.e. Hr = Hf = In the square bracket, the subtraction 
for Fff(x,fx) is defined as 



F c h (x^)sub = F 9 H {x,n)(£>F c g {x,n) = f F c g {n,n)F 9 H (-,?) —. 



(2) 



The quark distribution Fg(x,/i) inside an on-shell gluon up to order a s can be connected 
to the familiar g — > cc splitting function P g ^ c , i-e. Fg(x,/j,) = In -^Pg^^x), with 
'x) = \{1 — 2x + 2a: 2 ). Later on for convenience, we shall call the 'heavy quark mecha- 



g->c 



nisms', in which proper subtraction has been given according to method in GM-VFN scheme, 
as 'intrinsic ones' accordingly. 

In Eq.(l), the first term is the gluon-gluon fusion mechanism, the second and the third 
terms are the so called 'intrinsic' charm mechanisms [11], in which all the subtraction terms 
are necessary to avoid the double counting problem, since these terms represent the parts 
of the gluon-gluon fusion mechanism which are already included in a fully QCD evolved 
'intrinsic' charm distribution function [13]. The gluon-gluon fusion mechanism has been 
considered by several authors [4-8]. However, in these references, they usually used a PDF 
in a zero-mass variable-flavor-number scheme but performed the partonic cross section cal- 
culation using the non-zero heavy-quark masses. Such treatment shall not heavily affect the 



results for the hadronic production at LHC or TEVATRON as is the case of B c production 
[12], however it will make large discrepancies at the fixed target experiment, i.e. SELEX 
experiment. This is because, for the fixed target experiment, most of the generated H cc 
events are concentrated in the small p t regions, where large uncertainties are caused due 
to the inconsistent using of PDF. This is one of the reason that we adopt the GM-VFN 
scheme to study the hadronic production of H cc in which the heavy-quark mass effects can 
be treated in a consistent way both for the hard scattering amplitude and the PDFs. 

For the 'intrinsic' charm mechanisms at the leading order, we need to calculate subpro- 
cesses: g + c — > + c and c + c — > S cc . For the hadronic production via c + c — > H cc , 
because its hard subprocess is a 2 — > 1 subprocess and the p t cut is unavoidable to ensure 
the applicable of the PQCD calculation, it at least need to emit one hard gluon to obtain 
the p t distribution. Therefore we shall calculate c + c — > H cc + g other than c + c — > S cc in 
the following calculations. Note for c + c — > H cc + g mechanism, it has no double counting 
problem with the gluon-gluon fusion mechanism and does not need to introduce the sub- 
traction term, since it is one order higher than the gluon-gluon fusion mechanism according 
to the power counting rule shown in Ref.[13]. 

According to NRQCD formulation, the production rate of S cc can be factorized into two 
parts, one part is for the production of two or four free quarks (for the intrinsic mechanism 
or gluon-gluon fusion mechanism respectively) and is determined by pQCD, another part is 
for non-perturbative transition of the (cc)-diquark pair into S cc and can be defined in terms 
of non-relativistic QCD (NRQCD) [10] matrix elements. According to the discussions in 
Ref. [9], at the leading order of v c , the baryon S cc contains two configurations of the (cc)- 
diquark pair, one is that with the pair in (cc)g[ 3 yS'i], another is that in (cc) 6 [ 1 5'o]) whose 
matrix elements can be written as 

h x = ^(0|[^ ai e^ a2 +V a2 e^ ai ](a t a)^ a2t e^ ait |0), 
48 

h 3 = ^(0|[^ Ol e<rV Oa -V' a2 e^V Ol ](a t a)^ Oat ^eV' ait |0), (3) 

where aj(j = 1,2) label the color of the valence quark fields and a l {% = 1,2,3) are Pauli 
matrices, e = ia 2 . hi represents the probability for a (cc)-diquark pair in (cc) 6 [ 1 5'o] to 
transform into the baryon, while h 3 represents the probability for a (cc)-diquark pair in 
(cc)3[ 3 S'i] to transform into the baryon. According to the discussion in Ref. [9], both hi and 
h 3 are of order v 2 to |(O|xV'0| 3 -S'i)| 2 . The value of the two matrix elements hi and h 3 can 
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(a) 
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FIG. 1: The schematic Feynman diagrams for the hadroproduction of H cc from the gluon-gluon 
mechanism, where the dashed box stands for the hard interaction kernel. k\ and &2 are two 
momenta for the initial gluons, q c i and q C 4 are the momenta for the two outgoing c, P is the 
momentum of H cc . The (cc)-diquark pair is either in (cc)g [ 3 Si] or in (cc^^So] respectively. 

be determined with non-perturbative methods like QCD sum rule approach, however their 
values are unknown yet. The fragmentation of a diquark into a baryon is assumed to occur 
with unit probability and consequently, to have no influence on the production cross section. 
Further more, as the fragmentation function D(z) of a heavy diquark into a baryon, peaks 
near z ~ 1[6] 3 , and then the momentum of the final baryon may be considered roughly 
equal to the momentum of initial diquark. So to study the hadronic production of S cc is 
equivalent to study the hadronic production of (cc)-diquark. Under such condition, the 
value of NRQCD matrix element /13 can be naively related to the wave-function for the color 
anti-triplet [ 3 Si] cc state, i.e. = |\l/ cc (0)| 2 . And for convenience, since h\ and h s is of the 
same order of v c [9], we take hi to be /13 hereafter. 

The schematic Feynman diagrams for the gluon-gluon fusion mechanism are shown in 
Fig.(l). Fig.(l) shows that there are two ways for the two outgoing valence c quarks to 
form the (cc)-diquark pair and each way contains 36 Feynman diagrams that are similar to 
the case of hadronic production of B c (all the diagrams can be found in Ref. [18], and one 
only need to change all the b quark line there to the c quark line). However in Refs.[5- 
7], only Fig. (la) is considered and then only 36 Feynman diagrams have been taken into 
consideration. Since the contributions from the left and the right diagrams of Fig.(l) are 
the same and there is an f~) factor for the square of the amplitude by taking into account 
the symmetry of the diquark wavefunction, so there is an overall factor '2' for our total 

3 By taking a simple form of fragmentation function D(z), Ref.[7, 21] did a rough estimation for such effects. 
The results there show that such effect is really small. 
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(a) (b) (c) (d) (e) 




(f) (g) (h) (i) (J) 



FIG. 2: Typical Feynman diagrams for the sub-processes induced by 'intrinsic' charm. The upper 
five of them are those for g(jp\) + cfo) ^cc{pz) + c(p4) and the lower five are for c(pi) + c(p2) — ► 
^cc{pz)+g{P4) respectively. The (cc)-diquark pair is either in (cc)g [ 3 Si] or in (cc)e[ 1 So] respectively. 

cross-sections in comparing with those in Refs.[5-7]. In the present paper, as a cross check 
of the results in Refs. [5-7], we calculated it by using two different methods and made a 
cross check numerically between them. One method is to fully simplify the amplitude of the 
gluon-gluon fusion mechanism by using the improved helicity approach which was developed 
in case of the hadronic production of B c [18, 19]. More details of the calculation could be 
found in the appendix A. The other is to generate the Fortran program directly by the 
Feynman Diagram Calculation (FDC) program[20], which is a Reduce and Fortran package 
to perform Feynman diagram calculation automatically. The detailed treatment method of 
S cc in FDC could be found in appendix B. 

For the 'intrinsic' charm mechanism, we need to consider the following sub-processes, 
g + c — > E cc + c and c + c — > S cc + g, where the (cc)-diquark pair in S cc is in (cc)g[ 3 S'i] 
or (cc)e[ 1 5'o], respectively. Similar to the case of gluon-gluon fusion mechanism, There is a 
symmetry factor '2' for cross section. The typical Feynman diagrams are shown in Fig. (2). 
The final expression of the total square of amplitude is quit simple, and we adopt the FDC 
program [20] to obtain it directly. 

We will calculate the 'intrinsic' charm mechanism within the GM-VFN scheme. In GM- 
VFN scheme, when one talks about the heavy quark components of PDFs, and takes into 
account of both 'heavy quark mechanisms' and the gluon-gluon fusion mechanism for the 
hadronic production, one has to solve the double counting problem: i.e. a full QCD evolved 
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'heavy quark' charm/bottom distribution functions, according to the Altarelli-Parisi equa- 
tions, includes all the terms proportional to In {j^r^j (/-* the factorization scale and mg 
the heavy quark mass); and some of them come from the gluon-gluon fusion mechanism, 
i.e., a few terms appear from the integration of the phase-space for the gluon-gluon fusion 
mechanism. 

To be specific, according to Eq.(l), the inclusive S cc hadronic production via 'intrinsic 
charm mechanisms' can be formulated as, 

da = I ' dx i J dx 2 F t Hi (x 1 ,fi F ) x Fj l2 (x 2 ,fi)da ij ^ Scc x(x 1 ,x 2 ,fi) , (4) 

ij 

where % ^ j and i,j = g,c. Here, the heavy quark PDF F^(x,fi) (x = x 1 or x 2 , H = Hi 
or H 2 ), should include a proper subtraction term F^(x, h)sub as is defined in Eq.(2) in 
order to avoid the double counting of g + c — > S cc + g mechanism to the gluon-gluon fusion 
mechanism. da^^^ stands for the usual 2-to-2 differential cross-section, 

^ < 2X (27T) 4 |M| 2 A d* Pi (* \ 

da ij ^ ccX (x 1 ,x 2 ,^ ) = 2 2 11 (o^s(op\ S <L Pl ~ Pl ~ P2 ' ( 5 ) 

where p±, p 2 are the corresponding momenta for the initial two partons and ps, p^ are the 
momenta for the final ones respectively. The average over the initial parton's spins and 
colors and the sum over the initial and the final state's spins and colors are absorbed into 
\M\ 2 . All the expressions of |M| 2 , with all the mass effects being retained, are shown in the 
appendix B. 

The phase space integration can be manipulated by adopting the routines RAMBOS [22] 
and VEGAS [23], which is the same as that of the B c meson generator BCVEGPY [18, 19]. 



III. NUMERICAL CHECKS 



Before analyzing the properties for the hadronic production of S cc , we need to check the 
Tightness of program for all the mechanisms, especially, we should be more careful for the 
most complicate gluon-gluon fusion mechanism. 

First of all, all the programs are checked by examining the gauge invariance of the am- 
plitude, i.e. the amplitude vanishes when the polarization vector of an initial/final gluon is 
substituted by the momentum vector of this gluon 4 . Numerically, we find that the gauge 

4 All the Fortran codes are available from the authors on request. 
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FIG. 3: The energy dependence of the integrated partonic cross-section for the production of the 
baryons with heavy diquarks via the gluon-gluon fusion mechanism. The dotted line, solid line, 
dashed line and dash-dot line stand for the baryons with (cc)§ [ 3 <Si], (bc)^ [ Si], (bcj^So] and 
[ 3 5i] respectively. The curves for H cc and Ebb are divided by 2. 

invariance is guaranteed at the computer ability (double precision) for all these processes. 
Next, to make sure the rightness of our program for the gluon-gluon fusion mechanism, as 
mentioned before, the numerical results of our two programs agree with each other exactly. 

Furthermore, we compared our numerical results for the gluon-gluon fusion mechanism 
with those in the literature by using the same input parameters as were stated in the 
corresponding references. To make a complete comparison with the results listed in Ref.[5], 
we also calculate the partonic cross sections for the production of E bc and E b b through the 
subprocesses, gg — > Eb c + b + c with the (frc)-diquark in color-anti-triplet [ 3 Si] or I 1 So] state, 
and gg —>■ Sbb + b + b with the (fefe)-diquark in color-anti-triplet [ 3 5i] state. The programs for 
the production of E bc and S feb can be easily obtained from the program for the case of S cc . 

In Fig. (3), we show the partonic cross sections for the production of baryons with heavy 
diquarks via the gluon-gluon fusion subprocess. In drawing the curves, we adopt the same 
parameter values as were taken in Ref.[5], i.e. with a fixed value for a s (a s = 0.2) and 

|^ cc (0)| 2 = 0.039GeV 3 , |^> c (0)| 2 = 0.065GW 3 , |^ bb (0)| 2 = 0.152GeV 3 , (6) 

m c = 1.8GeV, m b = B.lGeV, M Hcc = 3.6GeV, M Sbc = 6.9GeV, M Sbb = !0.2GeV . (7) 

10 



TABLE I: Comparison of the partonic cross sections for gg — ► S cc + c + c with the corresponding 
results in Ref.[6], where the (cc)-diquark pair is in (cc)§ [ 3 Si]. E cm is the center of mass energy of 
the subprocess. The input parameters are m c = l.7GeV, M= cc = SAGeV, the radial wavefunction 
at the origin i? cc (0) = V4^^ cc (0) = 0.601Gey 3 / 2 and a s = 0.2. 





15 GeV 


20 GeV 


40 GeV 


60 GeV 


80 GeV 


100 GeV 


a{pb) 


66.6 


68.2 


41.8 


26.2 


17.9 


13.1 


a(pb)[6] 


23.2 


22.5 


13.7 


8.96 


6.45 


4.94 



For convenience of comparison with those of Ref.[5], in Fig. (3), our results for E cc and E bb 
have been divided by an overall factor '2'. One may easily find all the curves for the energy 
dependence of the partonic cross-sections shown in Fig. (3) are in consistent with the results 
in Ref.[5] (Fig. (2a) there). 

In Tab. (I), we show the comparison of partonic cross sections (the second column) for the 
production of E cc with the (cc)-diquark in (cc)g [ 3 Si] via the gluon-gluon fusion subprocess 
with those in Refs.[6, 7]. In Tab. (I), the results of Ref.[6] is derived from the fitted expression 
(Eq.(8) in Ref.[6]): 

° = 213. f 1 - V' 35 , (8) 



E \ E 

\ -'-'cm / \ -'-'cm / 

where E cm is the center of mass energy of the subprocess. One may observe that under the 
same parameter values, the results in Refs.[6, 7] are in disagreement with ours even though 
they are close in shape 5 . 

Next, as a cross-check between our results for one of the 'intrinsic' mechanism through 
c + c — > E cc + 9 with those of Ref.[21], we show the cross section of S cc -baryon production 
at the hadronic energy E cm = l.STeV or lATeV with the same input parameters in Fig.(4). 
The curves in Fig. (4) agree with those of Ref.[21] (Fig. 3 and Fig. 4 there). 

As a summary, for the gluon-gluon fusion mechanism, except for an overall factor '2', we 
confirm the results in Ref. [5], but not those of Ref.[6, 7]. And for one of the 'intrinsic' charm 
mechanism, i.e. c + c — > E cc + g, our results agree with those of Ref. [21] under the same 



5 Such discrepancy has already been found in Ref. [5] , however the author there attribute it to the different 
use of input parameters. 
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FIG. 4: The ^-distributions for the 'intrinsic' mechanism c + c — > H cc + g at Tevatron RUN-I 
E cm = 1.8TeV (dashed line) or LHC E cm = 14TeV (solid line) both with \y\ < 1.0. The present 
results are calculated with the same parameters in Ref. [21]. 
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FIG. 5: The Pt— and y— distributions of the produced E cc for the subprocess gg — > H cc + c+c under 
different C.M. energies. Here the (cc) diquark pair only in (cc)g[ 3 5'i] is taken into account. The 
dotted line, dashed line, solid line, dash-dot line and the diamond line stand for E cm = 20GeV, 
40GeV, 60GeV, 80GeV and lOOGeV, respectively. 



input parameters. 

Finally, we discuss the properties of the two different configurations of (cc)-diquark pair, 
i.e. (cc)^! 3 ^!] and (cc)^ 1 ^], for the hadronic production of S cc . In Fig. (5), we show the 
transverse momentum P t distribution and the rapidity y distribution at different center of 
mass energies for the subprocess gg — > E cc + c + c, with its (cc)-diquark pair in (cc^gpSi]. 
The case of (cc)-diquark pair in (cc^^So] is similar and will not be shown here. We drawn 
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FIG. 6: The energy dependence of the integrated partonic cross-section for the S cc production 
of via the gluon-gluon fusion mechanism. The solid line and the dashed line stand for the two 
(cc)-diquark pair states in configurations (cc)g[ 3 5i] and (cc^^So] respectively. 



a comparison between the energy dependence of the integrated partonic cross-section of the 
two different (cc)-diquark configurations, i.e. (ce)g[ 3 Si] and (cc) 6 [ 1 S , ], in Fig. (6). One may 
observe that the curves for (cc)g[ 3 Si] and (cc^^So] are close in shape and the contribu- 
tions from (cc)e [ 1 5'o] can be up to ~ 20% comparing with the case of (cc)g[ 3 »5'i]. So the 
contributions from the (cc)-diquark pair (cc^pSo] are significant and should be taken into 
consideration for a full estimation of the hadronic production of S cc . This is in agreement 
with the conclusion drawn in Ref. [9], where the contributions of these two different states 
of (cc)-diquark pair are discussed through the fragmentation approach. 



IV. HADRONIC PRODUCTION OF E cc 

In the present section, we shall first study the hadronic production properties of S cc both 
at TEVATRON and at LHC, and then make a discussion for the hadronic production at the 
fixed target SELEX experiment. All the calculations are done under the GM-VFN scheme. 
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TABLE II: Cross sections a for the hadronic production of H cc at TEVATRON and LHC, where 
the (cc)-diquark is in (ccjg [ 3 <Si] or (ccje^So], and the symbol g + c means g + c — > H cc + c and etc. 
In the calculations, p t > AGeV is taken and \y\ < 1.5 at LHC, while \y\ < 0.6 at TEVATRON. 





TEVATRON (y/S = 1.96 TeV) 


LHC (y/S = 14.0 TeV) 




(«0 S pSl] 




(cc)g[ 3 Si] 


(ccjet 1 ^] 


a g+g (nb) 


1.61 


0.392 


22.3 


5.44 


<j c+g {nb) 


2.29 


0.360 


22.1 


3.42 


a c+c (nb) 


0.751 


0.0431 


8.74 


0.475 



A. hadronic production of E cc at LHC and TEVATRON 

As has been discussed in Sec. II, we take h 3 = |\I/ CC (0)| 2 and hi = h 3 in the calculations. 
The mass of M= cc can be determined by potential model, and it is estimated to be[6], 
M Ecc = 3.584 ± 0.035GW. In Ref.[l], it has been measured to be 3.519 ± O.OOlGeV. For 
clarity, we choose |^ cc (0)| 2 = 0.039GeI/ 3 [5], M Scc = 3.50GeV and then m c = 1.75GeV. The 
factorization energy scale is fixed to be the transverse mass of S cc , i.e. Q — M t = \J M 2 + p 2 , 
where pt is the transverse momentum of the baryon. The PDFs of version CTEQ6HQ [16] 
and the leading order a s running above A Q ^ D = 0.32QGeV are adopted. 

In TABLE II, we show the cross-section for the hadronic production of S cc at TEVATRON 
and LHC, where p t > AGeV is taken in the calculations and \y\ < 1.5 at LHC, \y\ < 0.6 at 
TEVATRON. From Tab. II, one may observe that similar to the case of hadronic production 
of B c meson [12], the cross-sections of the 'intrinsic' charm mechanisms are comparable to, 
or even bigger than, the usual considered gluon-gluon fusion mechanism. From Tab. II, one 
may also observe that the contributions from (cc^^So] are sizable comparing with that 
of (cc)3[ 3 5'i], i.e. for the gluon-gluon fusion mechanism, the contribution from (cc)e P-S'o] is 
about 24% of that of (cc)^ 3 ^], while for the mechanisms of c+g — > H cc +c and c+c — > E cc +g, 
it changes to ~ 15% and ~ 5%, respectively. 

In Fig. (7), we show p t distributions for the hadronic production of S cc with two configu- 
rations of the (cc)-diquark pair states, i.e. (cc)g[ 3 Si] and (cc^^So], where \y\ < 1.5 at LHC 
and \y\ < 0.6 at TEVATRON are adopted. From Fig. (7), one may observe the following 
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FIG. 7: The ^-distribution for the hadroproduction of H cc at TEVATRON (left) and at LHC 
(right), where \y\ < 1.5 at LHC and \y\ < 0.6 at TEVATRON are adopted. The dotted line and 
the solid line are for gluon-gluon fusion mechanism, the triangle line and the diamond line are for 
g + c — * Sec, the dashed line and the dash-dot line are for c + c —* H cc , where the upper lines of 
each mechanism are for (cc^pSi] and the lower lines are for (cc^e^So], respectively. 

points: 1) to compare with the gluon-gluon fusion mechanism, the 'intrinsic' mechanism 
g + c —>■ S cc + g dominant in small p t regions and its p t distributions drop faster than that 
of gluon-gluon fusion mechanism, which is similar to the case of B c hadroproduction [12]. 
2) For 'intrinsic' mechanism c + c — > H cc + g, it pt-distribution drops faster than other 
mechanisms and then its contribution is the smallest among all the mechanisms. 3) For 
a particular mechanism, the contribution from the case of (cc)e P"*So] is sizable comparing 
with the contribution from the case of (cc)g[ 3 Si]. However, the p t distribution of (cc^pSo] 
is smaller than that of (cc)^ 3 ^] in the whole p t regions for the same mechanism and it 
also drops faster than the case of (cc)g[ 3 5'i]. Especially for the c + c — > E cc mechanism, 
Pt distribution of (cc^^So] drop much faster than that of (cc)3[ 3 5*i], and then the cross- 
section for (cc)e[ 1 5'o] is only about 5% of that of (cc)jj[ 3 5i]. As for the gluon-gluon fusion 
mechanism, the contribution from (cc^^So] is comparable to that of (cc)g[ 3 S'i] from the 
'intrinsic' mechanisms at high energies, especially at LHC, so one should be careful to take 
the contribution from (cc)e[ 1 >S'o] m t° consideration so as to provide a full estimation for all 
these hadronic mechanisms. 
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TABLE III: Cross section a for the hadronic production of H cc at the fixed target experiment with 
center of mass energy 33.58GeF, where the (cc)-diquark pair is in (ccJgpSi] or (cc^^So], and the 
symbol g + c means g + c — ► E cc + c and etc. In the calculations, pt > 0.2GeV is taken. 





SELEX (V£ = 33.58GeV) 








a c+c (pb) 


(cc)g[ 3 Si] 


4.03 


102. 


1.02 x 10~ 3 




0.754 


11.3 


4.15 x 10~ 5 



B. A simple discussion on hadronic production of E cc at the fixed target SELEX 
experiment 

For the fixed target experiment, the 'intrinsic' charm mechanism becomes more important 
than in the case of hadronic production at TEVATRON or LHC, since small p t events can 
contribute here. Such an experiment has been done by SELEX group [1] and it may cover 
all solid angle without p t cut, thus the 'intrinsic' charm mechanisms may be studied and 
extended to very small p t region. For SELEX experiment, its lower p t limit can be as small 
as 0.2GeV. However, one should be careful to ensure that the pQCD calculation is reliable in 
such small pt regions, i.e. the intermediate gluon (with momentum q) in all the mechanisms 
for the hadronic production of S cc must be hard enough, i.e. q 2 » A-qcd- 

For the gluon-gluon fusion subprocess, the square of the intermediate gluon momentum 
at least is bigger than (4m;?) so as to produce one cc-quark pair and then is always PQCD 
calculable. For the 'intrinsic' subprocess g{pi) + c(p 2 ) — > S cc (j9 3 ) + 0(^4), we must ensure 
that the momentum of the intermediate gluon of Fig. (2b, 2c, 2e) satisfy 

Q 2 = -q 2 = -(pi-^y»A 2 Q CD , (9) 
and for the 'intrinsic' subprocess c(pi) + c(p 2 ) — * S cc (p 3 ) + g(pi), similarly, we have 

Ql = -<ll = - (Pi - f ) 2 » Alco , Ql = -Ql = - (P2 - I)' » a qod- (io) 

Eqs.(9,10) give two extra constraints for both the partonic fractions x 1 , x 2 and p t . For 
definiteness, we set the lowest values for Q 2 , Q\ and Q\ to be ml (>> Aq CD ). 
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FIG. 8: The pt-distributions for the hadroproduction of S cc at SELEX. The dotted line and the 
solid line are for gluon-gluon fusion mechanism, the dashed line and the dash-dot line are for 
g + c — > H cc , the triangle line and the diamond line are for c + c — > H cc , where the upper lines of 
each mechanism are for (cc)jj[ 3 Si] and the lower lines are for (cc^e^So], respectively. 

We show the cross-section for the hadronic production of S cc at SELEX experiment in 
TABLE III, where pt > 0.2GeV is adopted in the calculations. TABLE III shows that at 
SELEX, the 'intrinsic' charm mechanism is the dominant mechanism and then the theoretical 
predictions of S cc events at SELEX can be raised by more than an order. We show the 
Pt distributions for the fixed target experiment in Fig. (8). One may observe that the pt 
distributions of 'intrinsic' mechanism g + c — > H cc + c are bigger than that of the gluon-gluon 
fusion mechanism almost in all the pt region, which is the reason why the total cross-section 
of g + c — > H cc + c mechanism is much larger than the gluon-gluon fusion mechanism as 
shown in TABLE III. For 'intrinsic' mechanism c + c — > H cc + g, it ^-distribution starts at 
~ 5GeV due to the constraint Eq.(lO) and its contribution is quite small. From TABLE 
III, one may also observe that the contributions from (cc^epSo] are also sizable comparing 
with that of (cc)g[ 3 Si] that is similar to the hadronic production at TEVATRON and LHC 
as shown in TABLE II, i.e. for the gluon-gluon fusion mechanism, the contribution from 
(cc^pSo] is about 19% of that of (cc)g[ 3 S'i], while for the processes of c + g — > S cc + c and 
c + c — > E cc + g, it changes to ~ 10% and ~ 4%, respectively. 
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TABLE IV: R values, which is defined in Eq.(ll), for the hadronic production of S, 





SELEX 
Pt > 0.2GeV 


TEVATRON 

Pt > 4GeV, \y\ < 0.6 


LHC 

Pt > 4:GeV, \y\ < 1.5 


R 


29. 


3.4 


2.8 



V. SUMMARY 

We have calculated the hadronic production of the doubly charmed baryon S cc via 
the gluon-gluon fusion mechanism and the 'intrinsic' charm mechanism, i.e. via the sub- 
processes g + g — > S cc + • • • , g + c — > S cc + • • • and c + c — > S cc + • • • . To avoid the double 
counting problem while taking the gluon-gluon fusion mechanism and the 'intrinsic' charm 
mechanism into consideration, we have adopted the GM-VFN scheme in which the heavy- 
quark mass effects can be treated in a consistent way both for the hard scattering amplitude 
and the PDFs. Some cross checks for the present results with those in the literature have 
been done. The result for the gluon-gluon fusion mechanism agree with what was given in 
Ref. [5] up to a factor of 2; and the results for the c + c — > S cc + g with (cc)-diquark pair in 
ccgpS'i] agree with that of Ref. [21] when adopting the same input parameters. Whereas 
the results for the 'intrinsic' mechanisms and those for the cases with (cc)-diquark pair in 
cc^So] are fresh. 

From TABLE II and TABLE III, one may see that the total cross sections of the 'intrinsic' 
charm mechanisms are comparable to, or even bigger than, that of the gluon-gluon fusion 
process, especially for the g + c — > S cc mechanism. To be more definite, we define a ratio 

R = ^ , (11) 

0- flfl ^S cc ((cc)3[35i]) 

where a tota i stands for the cross section for all the concerned mechanisms and o" 9 (,-»s cc (( cc )3[ a 5i]) 
is the cross section for the gluon-gluon fusion mechanism with (cc)-diquark pair in (cc)g[ 3 S'i] 
configuration only. The values of R for the hadronic production of S cc in various environ- 
ments are shown in Tab. IV, which shows that the 'intrinsic' charm mechanisms are not 
negligible: at SELEX they even dominate over the other mechanisms. The contributions 
from the (cc)-diquark pair in cc^So] for all the concerned mechanisms are also considered 
in the work, and the results show that if the matrix element hi is at the same order of h 3 [9] , 
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i.e. hi ~ h,3, the diquark pair will make a sizable contribution to the hadronic production 
of E cc . 

We may conclude that to be a full estimation for the hadronic production of H cc , one 
needs to take all these mechanisms into consideration. One may observe that by taking 
into account the 'intrinsic' mechanisms, the theoretical prediction on the S cc event can be 
almost one order higher than the previous predictions in which only the gluon-gluon fusion 
mechanism is considered. Nevertheless, there is still a big discrepancy between the SELEX 
observation [1] and pQCD predictions. Perhaps it is due to the fact that the small p t regions 
is not amenable to the pQCD analysis, e.g., the 'intrinsic' mechanism c + c — > (cc)[ 3 Si]3 + g 
and c + c — > {cc)[ 1 Sq]q + g, according to constraint (10) there is a big contribution from non- 
perturbative QCD range, and the intrinsic charm fusion mechanism with the subprocesses 
c + c — > (ccJpSijg and c + c — > (cc) p-Sb] 6 j which may contribute to the production greatly 
but only with very small p t , and being non-perturbative QCD nature, it is not considered 
here. Another possibility might be that the SELEX experiment does not provide sufficient 
support for its claim of evidence for the observation of doubly charmed baryon S cc as pointed 
out by Ref.[3]. 

Acknowledgments '. This work was supported in part by the Natural Science 
Foundation of China (NSFC). 

APPENDIX A: CALCULATION TECHNOLOGY FOR THE GLUON-GLUON 
FUSION MECHANISM UNDER THE IMPROVED HELICITY APPROACH 

The general structure of the amplitude in 'explicit helicity' form can be written as 

M^ 2 ' XM \q c3 , g c4 , q cl , q c2) k 1: k 2 ) = g A s £ W^iJ^^'^'fe, g c4 , q cl , q c2 , h, k 2 ) ■ 

A2,A3 

D 2 B^\ qc3 ,q cl ), (Al) 

where % — 1, • • • , 72, Aj (j — 1, • • • , 6) denote the helicities of the quarks and gluons respec- 
tively. Ai denotes the helicity of c(q C 3), X 2 that of c(q C 4), A3 that of c(q c i), A4 that of c(q c2 ); 
whereas A 5 denotes that of gluon-1 and A 6 denotes that of gluon-2. Here Cj, Xi denote the 
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color factor and the scalar factor from all the propagators as a whole for the ith-diagram, 
respectively. B^l MMMM ' X6 \q c3 , q c4 , q cl , q c2 , h, k 2 ) and B^ Xs \q c3 , q cl ) are the amplitudes 
corresponding to the 'free quark part' g(ki, \§)g{k 2) A 6 ) — > c(g c3 , Ai) + c(g c4 , A 2 ) +c(q c i, A 3 ) + 
c(q C 2, A4) (all the quarks are on-shell) and the 'bound state part' c(g c3 , Xi) + c(q cl , A 3 ) — > (cc), 
respectively. L>i = /r) 1 ,„ 1 /0 1 /0 1 and _D 2 = m 1 /o 1 are two common 
normalization factors. 

By comparing Eq.(Al) with Eq.(22) in Ref.[18] that is for the B c hadropro duct ion, one 
may observe that both amplitudes are quite similar with each other. Most of the present 
helicity amplitudes can be directly derived from the results in Ref. [18] by simply replacing 
the 6-quark line there to the present c-quark line. And for the present case, we only need 
to deal with the following type of the helicity matrix element (HME) that is quite different 
from the case of B c hadroproduction, i.e. 

HMEj = {qoxMa + m c )t ^ - m c )\q 0Xl ) , (A2) 

where % = (1, • • • , 72) stands for the i-th Feynman diagram and f ; means that all the 
momentum in Tj (Tj stands for the explicit strings of Dirac 7 matrices between U (g c3 ) and 
V(q c4 ), which corresponds to i-th Feynman diagram) should change their sign and the string 
of the 7-matrices in Tj should be written in inverse order. In fact, such type of HME can 
also be relate to the familiar one as has been dealt with in the B c case by adopting the 
following relation: 

HMEj = -(<7o(-Ai)Wc3 + m c )Yi(i cA - m c )|g 0( _ A2) ) . (A3) 

A simple demonstration of Eq.(A3) can be found in the last part of the appendix. 

The sum of all the helicity amplitudes of the sub-process g + g — > (cc) + c + c can be 
arranged as 

6 

M^ XM *\q c3 ,q cA ,q cl , q c2 ,h,k 2 ) = ^ C rnij M^ 2 ' > " i ' X6 ' X6 \q c3 , q c4 , q cl , q c2 , ki, k 2 ) ,(A4) 

m=l 

where C mi j (m — 1 — 6) are six independent color factors of the process, 

(~i ' (rparpb\ f~i /~1 ' (rpbrpa\ f~i 

I I 

C 3 ij 2^\/2 ^ m ^^ ] )niG m nk i C 4 ij 2,\^2,^^ ^ Tra ^^ _ ' ) n *^ mn ^ ' 
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TABLE V: The square of the six independent color factors (including the cross terms) for gg 
(cc)g[ 3 5i] + c + c, (C mi j x C*j ■) with m, n = (1, 2, ■ • • , 6), respectively. 
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TABLE VI: The square of the six independent color factors (including the cross terms) for gg — ► 
(ccJepSo] + c + c, (Cmi,- x C*j ■) with m, ra = (1, 2, • ■ ■ , 6), respectively. 
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where i,j = 1,2,3 are color indices of the two outgoing anti-quarks c and c respectively, 
and the indices a and b are color indices for gluon-1 and gluon-2 respectively. Here, the 
function G m jk equals to the anti-symmetric e m jk for the (cc)-diquark in 3 configuration and 
equals to the symmetric f m jk for the (cc)-diquark in 6 configuration respectively. The anti- 
symmetric e m jk satisfies e m jk^m'j'k = Smm'Sjj' — bmj'bjm' an d the symmetric f m jk satisfies 

fmjkfm'j'k ^mm'^jj' Omj'vjm' • 

To get the matrix element squared, one needs to deal with the square of the above six 
independent color factors as shown in Eq.(A5) (including the cross terms), i.e. (C mi j x C*^) 
with m,n = (1, 2, • • • 6). For reference use, we list the square of these six independent color 
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factors in TABLE V and TABLE VI, which are for (cc)g[ 3 S'i] and (cc) 6 [ 1 S , ], respectively. 

By keeping all these points in mind, we rewrite a program based on the B c meson gen- 
erator BCVEGPY[18, 19] to calculate the gluon-gluon fusion mechanism for the hadronic 
production of S cc . 

Finally, we give a simple demonstration of the relation Eq.(A3). To demonstrate the 
relation Eq.(A3), we shall adopt the following relation, 

(P(M)\h-h\q ( x 2 )) = {-l) n+1 {q^x 2) \h-4M-x,)) , (A6) 
whose non-zero ones can be explicitly written as [25] 

(p_|#i...# n |g + ) = -(q-\lt n ...jt 1 \p + ) (n even), (A7) 
(p+|#i...# n |g_) = -(q+\h~-h\P-) ( n even), (A8) 
(p+\h-h\q+) = (q-\h-h\P-) in odd), (A9) 

where ki(i — 1, • • • , n) are any types of momenta. 

Generally, to the i-th Feynman diagram, we can expand Tj as, 

n 

and then we have, 

ft = £(-i) B <W„--- AA), 

n 

where C n are functions free of Dirac 7 matrix element. Taking use of Eq.(A6), 
obtain 

(<7o(-Ai)Wc3 + m c )T 2i (<icA - m c )|g (-A2)) 

n 

+ m c (^ ^ 2 • • • jf n ) 4 cA - m 2 c {^ i> 2 --- i>^\ |go(-Aa)) 

= £ Cn(?0(A 2 )| [(-1)" +3 (U^ • • ' A A) 4c3 ~ (-l)" +2 m c (^ • • • i> 2 A) 4c3 

n 

+ (-l)" +2 m c ^ •••^A)- ■ • • f 2 A)] \q 0{Xl) ) 

= -{qo(\2)\(4c4 + m c )f2i(4c3-m c )\qo( Xl ))- (A12) 
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APPENDIX B: CALCULATION TECHNOLOGY IN FDC PROGRAM[20] AND 
THE SQUARE OF AMPLITUDE FOR THE INTRINSIC CHARM MECHANISM 



First, we take gluon-gluon fusion mechanism as an explicit example to show the technol- 
ogy used in FDC program[20] and show in more detail how we can derive the program for 
the hadronic production of S cc from those of J/ip. 

The amplitude for each Feynman diagram of g + g — > J/ipipz) + c(p 4 ) + c{p<j) can be 
written as: 

- P3 

M(J/ip) = u(p^s 4 )T l s f (k l ,m c ) ■ ■ ■ S/(fc n _i, m c )T n v( — , si) 

P3 

B(p 3 , s, si, s 2 , mj/^)u(—, s 2 )T' 1 s f (q 1 , m c ) ■ ■ ■ s/(g„/_i)r / n ,u(p 5 , s 5 ). (Bl) 

where Sf(k,m) (k = fc, or q{) is the fermion propagator, B(p 3 , s, s±, s 2 , ttlj/^) is the wave- 
function of J/ip, Ti, • • • , r„, r' 1; • • • , F' n , are the interaction vertices. The color factor part is 
treated separately (similar to the method described in Appendix.A) and will not discussed 
here. 

One can easily find out the corresponding Feynman diagram in g + g — > S cc (p 3 ) + c(p 4 ) + 
c(p 5 ) and the amplitude of it could be written as: 

P3 

M(E CC ) = u( — ,si)T n s f (-k n -i,m c )---Sf(-k u rn c )T 1 v{p 4 ,s A ) 

P3 

B(p3,s,s u s 2 ,m^ cc )u(—,s 2 )Y' l s f (q 1 ,m c ) • • • s/(<V-i, m c )V n ,v(p 5 , s 5 ), (B2) 
where B(ps, s, s±, s 2 , ^s cc ) is the wavefunction of S cc . For an arbitrary Fermion line, 

P3 

a = u( — , si)r n s/(-A; n _i, m c ) ■ ■ ■ s/(-/ci, m^T^v (p 4 , s 4 ), 

we have 

a = a T = v T (p 4 , s 4 )rf sj(-/ci, m c ) • • • sj(-fc n _i, m c )r^«(y , si) 

= w T (p 4 , S4)CC , -r[CC-sT(-A; 1 , m c )CC- • • • CC-sT(-jk n _i, m c )C , C , -rJTC'C'-u(^, Sl f 

2 

= (-l) (n+1) M(p 4 , S 4 )riS/(fci, m c ) ■ ■ ■ Sf(kn-!, m c )r n w(y , Si), 

with the help of the following equations 

V T (p 4 ,S 4 )C = -U(p 4 , S 4 ), C~u(y,Si) =u(y,Si), 

CTf C = -r,, C-s^-fci, m c )C = s / (fc i , m c ). 
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Where C = —i^y 2 ^ is the charge conjugation matrix. And then Eq.(B2) can be transformed 

as 

M(S CC ) = (-l) (n+1) w(j»4, s 4 )T 1 s f (k l ,m c ) ■ ■ •s / (fc n _i,m c )r n t;(y,si) 

P3 

B(p 3 , s, si, s 2 , m Ecc )u( — , s 2 )T , 1 s f (q 1 ,m c ) ■ ■ ■ s f {q nl - 1 )Y' n> v{p b , s 5 ). (B3) 

By comparing Eq.(Bl) with Eq.(B3), one find that they are the same except for an overall 
factor (— l)( ra+1 ) ? where 'n' is the interaction vertex number of the corresponding fermion 
line and depends on the detail of each Feynman diagram. Therefore, we can completely use 
the method of J/ip to deal with S cc case by adding a factor (— l)( n+1 ) diagram by diagram. 
The detailed description of method to treat the J /if; and B c calculation could be found in 
the Ref.[19, 20]. 

All the above discussion is also valid for the calculation of the 'intrinsic' charm mecha- 
nisms. And the following results are obtained by taking the FDC program. 

For convenience, we express the square of the amplitudes by the variants s, t and u, which 
are defined as: 

s = (pi + p 2 f, t = (pi - p 3 ) 2 , u = (p!- p 4 ) 2 , 

where pi = (Ei,p ix ,pi y ,pi Z ) are the corresponding momenta for the involved particles: p\ 
and P2 are the momenta of initial partons, p 3 and p^ are the momenta of H cc and another 
outgoing particles respectively. Further more, for c(pi) + g{p 2 ) —> ^ccipz) + c{p±), we set 

ui — (u — Am 2 c ), s 1 — (s — m 2 c ), ti — (t — ml), 

and for c(pi) + c(p 2 ) -> S cc (p 3 ) + g(p 4 ), we set 

Ui — (u — ml), si — (s — Am 2 c ), h — (t — m 2 c ). 

The relation, u\ + 1\ + s\ = 0, is useful to make all the expressions for the square of the 
amplitudes compact. 

The square of the amplitude for the subprocess c(pi) + g{p 2 ) — > 2 cc (p 3 ) + c{p±) with 
(cc)-diquark pair in (cc)g [ 3 Si] can be written as, 
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The square of the amplitude for the subprocess c(pi) + g{j>2) — ► 2 cc (p 3 ) + c(p 4 ) with 
(cc)-diquark pair in (cc)e[ 1 5'o] can be written as, 



m = 2_W|^ cc (0)| 2 7r 4 
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The square of the amplitude for the subprocess c(pi) + c(p 2 ) 
(cc)-diquark pair in (cc)s [ 3 Si] can be written as, 
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The square of the amplitude for the subprocess c(p\) + c(p2) 
(cc)-diquark pair in (cc) 6 [ 1 S'o] can be written as, 
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In these equation, M is the mass of H cc and \I/ CC (0) is the wavefunction at origin for the [ 3 Si] 
cc state. And here we have adopted that h 3 = |\l/ cc (0)| 2 and hi = h 3 . 
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